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ABSTRACT 

We present the latest velocities for ten multiplanet systems, including a re- 
analysis of archival Keck and Lick data, resulting in improved velocities that 
supersede our previously published measurements. We derive updated orbital fits 
for ten Lick and Keck systems, including two systems (HD 11964, HD 183263) 
for which we provide confirmation of second planets only tentatively identified 
elsewhere, and two others (HD 187123 and HD 217107) for which we provide 
a major revision of the outer planet's orbit. We compile orbital elements from 
the literature to generate a catalog of the 28 published multiple-planet systems 
around stars within 200 pc. From this catalog we find several intriguing patterns 
emerging: 



• Including those systems with long-term radial velocity trends, at least 28% 
of known planetary systems appear to contain multiple planets; 

• Planets in multiple-planet systems have somewhat smaller eccentricities 
than single planets; and 

• The distribution of orbital distances of planets in multiplanet systems and 
single planets are inconsistent: single- planet systems show a pileup at P ~ 3 
days and a jump near 1 AU, while multiplanet systems show a more uniform 
distribution in log-period. 

In addition, among all planetary systems we find the following: 

• There may be an emerging, positive correlation between stellar mass and 
giant-planet semimajor axis. 

• Exoplanets more massive than Jupiter have eccentricities broadly dis- 
tributed across < e < 0.5, while lower mass exoplanets exhibit a dis- 
tribution peaked near e = 0. 

Subject headings: planetary systems 



Introduction 



1.1. The Detection of Multiple-Planet Systems 

Th e first exoplanetary system known to comprise multiple planet^ was the triple system 
V And (jButler et al.lll999l ). detected by the radial velocity m ethod just four years after the 
first confirmed exoplanet, 51 Peg b ( Mavor fc Quelo3 19951 ). The subsequent discovery of 



a second planet orbiting 47 UMa (Fischer et al.l l2002al ) and the resonant pair of planets 



orbiting G J 876 (IMarcy et al.l l200ll ) foreshadowed the discovery of more than two dozen 
more systems. Today, 14% of known host stars of exoplanets within 200 pc are known to 
be multiple-planet systems, and another 14% show significant evidence of multiplicity in the 
form of long-term radial velocity trends. 

The recent proliferation of multiple-planet systems is due to the increase in both the 
velocity precision and duration of the major planet search programs. The increased time 
baseline has led to the detection of long-period outer companions, the first being the P = 



^Prior to this, Wolszczan fc Fraill (1992) detected three extraordinary planets orbiting the pulsar PSR 
1257+12. Here, we restrict the discussion to systems orbiting nearby, ordinary stars. 



1270 d planet in the multiple system v And. Today, the 10+ year baseline of high-precision 
(< 5 m/s) radial velocity planet searches means that most planets with a < 3 AU now 
have multiple complete orbits observed, improving their detect ability. Even planets with 
a > 4 AU which have not yet completed a single orbit can sometimes have we ll constrained 
minimum masses (Msini), as in the case of HD 187123 c (jWright et al.l 120071 ). 



Radial velocity precision has steadily improved toward and below 1 m/s with the High 
Resolution Echelle Spectr ometer (HIRES) at the Keck Observ atory and with the HARPS 
spectrograph at La Silla (jPepe et al.l l2003l : IWright et al.l 120071 ). This has allowed for the 
detection of ever weaker signals and led to t he discovery of some of the lowest- mass planets 
known, including (among others): 55 Cnc e (McArthur et al.l 12004 ) , GJ 876d (IRivera et al. 



20051 ). the triple system HD 69830 (ILovis et al.ll2006l ). and /x Ara d (ISantos et al.l 12004 ). We 
list all published multiple-planet systems within 200pc in Table [H a list which benefits from 
updated and improved velocities of 10 of the systems provided in Table [2l We present a 
graphical overview of these multiplanet systems in Fig. [H 



Table 1. List of Exoplanets in Multiplanet Systems 
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Fig. 1. — Chart of semimajor axes and minimum masses for the 28 known multiplanet 
systems. The diameters depicted for planets are proportional the cube root of the planetary 
Msini. The periapse to apoapse excursion is shown by a horizontal line intersecting the 
planet. The diameters depicted for stars are proportional the cube root of the stellar mass. 
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Note. — For succinctness, we express uncertainties using parenthetical notation, where the least significant digit of the uncertainty, 
in parentheses, and that of the quantity are understood to have the same place value. Thus, "0.100(20)" indicates "0.100 ± 0.020", 
"1.0(2.0)" indicates "1.0±2.0", and "1(20)" indicates "1±20". 



^Whcn the uncertainty in e is comparable to e, uncertainties in uj and e become non-Gaussian. See lButler et al.l (|2006r ) for details. 

Bean et all (Eooj); Bu 6: iButler et al 



^References indicate which o rbital paramet ers are taken f r om th e literat ure as follows: Be8: 



(|2006f ): Cc7: ICochran et all (|2007l ): Cr5: ICorreia et al.l (|2005l ): Ds8: iDesort et al.l (|2008l ): Fi2: iFischer et all (|2002b[ ): Fi8: 



Fischer et alJ 120081): L e6: iLee et all (I2OO6I): Ly 6: ILovIs et all (I2OO6I): M v4: iMavor et all (I2OO4I): Mv8: Mayor et all (I2OO8I) : Ni8: 
Niedzielsk i et al.' ('2OO8I): Nf4: iNaef et all (I2OO4I); Pp7: IPepe et all ( 20071): R5: iRivera et all |2005h : T6: Ixinnev et all |2006l) : U7: 
Udrv et al. (2007 ): Vo5: IVoet et al.l (|2005l ): Wr7: IWright et al.l (|2007l)~ All other orbital solutions are new Keplerian (kinematic) fits 
to the data in Table [H 

'^Tho fit for this system includes a trend of —0.51 ± 0.1 m s^^yr^^. 

'^The planets in HD 73526 are in a 2:1 mean motion resonance, and planet-pl a net in teractions arc important, rendering Keplerian 
elements inadequate. In addition to the elements reported here, Tinnev et al. ( 2006h report a mean anomaly to be 86° ± 13 and 
82° ± 27 at a Julian Date of 2451212.1302. There is considerable degeneracy between K and e because the orbital period of HD 
73526c differs from one year by only 12 d. 

^ Barnes. Gozdziewski. fc: Ravmondl ( 20081 ) found that the orbit presented in iBean et al.l ( 2008h for the d component is unstable. 



and provide multiple stable solutions without uncertainties in the orbital parameters. 

^Planet-plan et interactions are st rong in 55 Cnc. The osculating orbital elements here are from the dynamical fit at Julian Date 



2447578.730 of Fischer et al 



( 2008h . That work puts no errors on these parameters, however, so the errors quoted here are those 
from the Keplerian (kinematic) fit there. 

^The exoplanets in HD 82946 have significan t interacti o ns, wh ich renders Keplerian orbital elements inadequate for describing their 



Lee et al 



(200 



B) report the mean anomaly of the inner and outer planets to be 353° 



orbits, since these elements are time-variable, 
and 207°, respectively, at a Julian Date of 2451185.1. 

'^This solution includes a linear trend with magnitude —3.08 ±0.16 m s~^ yr~^. 

'The exoplanets in HD 202206 have significa nt interactions , whic h renders Keplerian orbital elements inadequate for describing their 
orbits, since these elements are time- variable. ICorreia et al.l (j2005| ) report the mean longitude to be 266.23° ± 0.06 and 30.59° ± 2.84 
for the inner and outer planets, respectively, at a Julian Date of 2452250. 

J The outer two exoplanets GJ 876 have significant in teractions, wh i ch ren ders Keplerian orbital elements inadequate for describing 
their orbits, since these elements are time-variable. iRivera et al.l (|2005l ) report the mean anomaly of the planets to be Md = 
309.5° ± 5.1°, Afc = 308.5° ± 1.4°, and = 175.5° ± 6.0° respectively, at a Juhan Date of 2452490. The solution quoted here 
assumes i = 90°. 



1.2. Multiple-Planet Systems and Planet Formation Theory 



multiple-planet systems are of special interest to test theoretical models of planet for- 
mation, dynamics, and final architectures. To date, over 230 planets have been discovered 
orbiting 205 main sequence stars within 200 pc, 14% of which harbor multiple planets with 
well constrained masses and periods. With 67 such planets members of multiplanet systems, 
we may now make statistically significant comparisons between properties of planets in single 
planet systems with those in multiplanet systems. 

The known planet population is remarkably diverse with properties that bear on planet 
formation theory. Theories to explain the semimajor axis distributions and eccentricity dis- 
tributions of planets have especially benefitted from the constraints that multiplanet systems 
impose. Most planet formation theories are based on the core accretion model that begins 
with a disk of dust and gas where the du s t particles collide and grow to form rock-ice plane- 
tary cores (lAarseth. Lin, fc Palmer! 1 19931 : iKokubo fc Idall2002l : iLevison. Lissauer. fc Duncan 
19981 ). If a core becomes massive enough while gas r emains in the disk, it gravitationally 



accretes the nearby gas and rapidly increases in mass (IBodenheimer. Laughlin. fc Liru 12003 



Pollack et al.lll996l ). Such gas planets should form preferentially beyond the "ice line" (near 
3 AU for solar- type stars), where ices can participate in the initial planetary cores. This 
appears inconsistent with the observation that about 20% of known exoplanets orbit within 
0.1 AU, where there should be too little ice in the protoplanetary disk for massive cores to 
form quickly. 

Thus it appears that short-period planets form farther out and migrate inwards to 
their final semimajor axis (jPapaloizoul l2005l : iTanaka fc Wardll2004l ; iTrilling. Lunine. fc Benz 
2OO2I ). The discovery of at least five systems in or near mean-motion resonance (viz. GJ 876, 
55 C nc, HP 82943, HP 73526, and HP 128311) may lend support to the migration hypoth- 
esis ( jNelson fc Papaloizoul l2002l : iLee et al.l l2006l : iMarzari. SchoU. fc Tricaricd l2005l ) . Such 
resonances are difficult to explain if planets form in situ, b ut hydrodynamical simulations 



and n-body simulations with externally applied damping (jKley et al. 



2005 



Brvden et al. 



2OO0I : IChiang fc Murravl I2OO2I : IP'Angelo Kiev, fc Hemiin3l2003l : llda fc Linll2004[ ) show that 



resonance capture occurs if planets undergo significant migration at different rates, passing 
through mean-motion resonances. Also, although there are confirmed 2:1 mean motion res- 
onances, other period ratios may occur (IKley. Peitz. fc Bryderu |2004J : iNelson fc Papaloizou 
2OO2I : iLaughlin &: Chambers! l2002l ). Future discoveries of multiplanet systems may shed fur- 
ther light on the resonance capture processes and on planet formation dynamics. 

multiplanet systems also provide hints about the wide distribution of exoplanet ec- 
centricities and we are beginning to identify different processes that drive so me planets to 



large eccentricities, while damping others to moderate and low eccentricities ( iMurrayl 12003 



Table 2. Updated RV Data for Multiple- Planet Systems 



Star Time Radial Velocity Unc. Telescope 







(JD-2440000) 


(m/s) 


(m/s) 
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And 


9680.753969 
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6.8 
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upsilon 


And 


9942.007812 
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9.1 
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upsilon 


And 


9969.979687 


-22.4 


6.5 
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upsilon 


And 


9984.852758 


-55.0 


6.9 
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upsilon 


And 


10032.741397 


28.7 


8.0 
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upsilon 


And 


10056.842667 


47.9 


6.4 
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upsilon 


And 


10068.585579 


-60.2 


7.1 
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upsilon 


And 


10068.773435 


-62.9 


7.4 
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upsilon 


And 


10069.598633 


2 


12 


L 


upsilon 


And 


10072.602539 


-91 


10 


L 



Note. — A full version of this table is available in the electronic 
version of the Journal This is only a sample. 



Barnes fc Quinru |2004| ). iFordI (120061 ) suggests that planet-planet scattering augments ec- 
centricities, and that after the era of strong planet-planet scattering interactions with the 
remaining planetesimals will damp these eccentricities to the observed distribution. One may 
hope to detect some signature of these process in exoplanet data. Indeed, the configuration 
of the V Andromedae three-planet system may be explained by a planet-planet scattering 
event that ejected a fourth planet from the vicin ity, as this system carries a signature of 
sudden perturbation (jFord. Lystad. fc Rasid l2005l) . HP 128311, a two -planet system, also 
appears to carry a signature of similar scattering (ISandor fc Kleyll2006l ). 



Models of planet formation in protoplanetary disks point to the disk viscosities and 
lifetimes that are required to produce some observed systems. There is some inconsis- 
tency between the required d isk behavior and that predicted by hydrodynamical theories 
( iKley. Peitz. fc Bryderu l2004l ). Simulations that include multiple protoplanets interacting 
with a protoplanet swarm show that most protoplanets will ev entually accrete on to the star, 
leading to a planet occurrence rate lower than that observed (ICresswell fc NelsonI 120061 ). 



The increasing number of characterized multiplanet systems should help constrain these 
theories. We present the current catalog of multiple-planet systems with the most recent 
radial velocity data in hopes of offering theorists a thorough overview current state of obser- 
vations, and an opportunity to propose new tests of planet formation theories as the number 
of known multiple-planet systems grows and our knowledge of them improves. 



1.3. Plan 

In §2] we present the sample of the 28 published nearby multiple-planet systems, includ- 
ing updated orbital parameters and radial velocities for ten of these systems. In [|3]we present 
some empirical correlations among the properties of multiple-planet systems and difference 
between the distributions of orbital parameters among multiple-planet and single-planet 
systems. In ||l]we briefly summarize our findings and discuss our principle conclusions. 



2. A Survey of the Known Multiple- Planet Systems 

2.1. New and updated multiplanet systems 

We present updated velocities and orbital parameters for ten multiple-exoplanet systems 
in Tables [Hand [21 In six systems, the updated orbits are substantially similar to the published 
orbits. In four cases there are major updates to previously published orbital parameters. 



which we detail below. 



2.1.1. HD 183263 



Marcy et al.l (120051 ) announced an Msini = 3.7Mj^p planet in a 634-day, eccentric 
0.4) orbit around HP 183263, a G2IV star, and pointed out a strong residual linear 



trend of 32 m s ^yr ^. IWright et al.l (120071 ) showed that by 2007 the residuals had signifi- 



cant curvature, but could not constrain the minimum mass of the outer companion. Since 
then, we have obtained observations on an additional six nights. These data have dramat- 
ically constrained the orbital fit to the outer planet, and we can now confidently measure 
the orbital elements of a 3-4 Mj^p planet under the assumption that t here is not a third 
companion contributing a detectable linear trend (see lWright et al.ll2007l . for a discussion of 
the difficulties of constraining planets with only partially-observed orbits). HD 183263c has 
a 8.4±0.3 yr orbital period, orbits at 4.3±0.4 AU, and has e = 0.24 ± 0.06. We present the 
full set of newly determined orbital parameters in Table [T], and the latest radial velocities in 
Table [2] and Fig. [2j Figure [3] shows in P — Msini space, demonstrating that, despite our 
having only one observed orbit, the minimum mass and orbital period of the planet are well 
constrained. 

In Table [1] we report a stable solution near the minimum with errors est i mated from 
the sample of stable orbits found through error bootstrapping (IWright et al.l 120071 ). We 
tested each of the 100 bootstrapping trials for stability, and found a rough dividing line 
between stable and unstable solutions for this system such that orbital solutions with Cc > 
10~^((Pc/days) — 3000) are generally unstable. We have mapped this line into Pc — irtcsaiic 
space in Figure [3] demonstrating that many solutions consistent with the data are, in fact, 
unstable. The true uncertainties are thus asymmetric about the nominal values because the 
minimum is so close to the boundary of stability. 

For these long-term stability tests, we applied direct n-body integrations on each of 
the orbital solutions generated in the error bootstrapping. The radial velocity parameter s 
were converted into initial conditions using a Jacobi coordinate system (ILee fc Peald l2003l ) . 
Unless otherwise specified, ^ ye assumed edge-on, coplanar orbits. We held the stellar mass 
fixed, adopting values from iTakeda et al.l (|2007n . We integrated for at least 10^ yr using 
the hybrid integrator in Mercury ( ChambersI 1999 ). For the majority of each integration. 
Mercury uses a mixed- variable symplectic integrator (| Wisdom &: Holmanlll99ll ) with a time 
step equal to a hundredth of the Keplerian orbital period calculated at a semimajor axis 
equal to the pericenter distance of the closest planet. During close encounters. Mercury uses 
a Bulrich-Stoer integrator with an accuracy parameter of 10"^''. We identified each set of 
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Fig. 2. — RV curves for HD 183263. The data are from Keck Observatory and show (top) 
the inner planet with P — 1.7 yr and (bottom) Msini = 3.7Mj^p the outer planet with 
P ~ 8 yr and Msini = 3.6 Mj^p. 
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Fig. 3. — Contours of in — rric sin i space of best-fit orbits to the RV data of HD 183263 
(Fig. [2]), with in grayscale. The sohd contours mark the levels where increases by 1, 4 
and 9 from the minimum. The dashed curve marks the approximate dividing line between 
stable and unstable orbital solutions, as determined by an ensemble o f n-body integ r ations of 
the system assuming edge-on, coplanar orientations (cf. Figure 13 of IWright et al.l (120071 )). 



initial conditions as an unstable system if: 1) two planets collide, 2) a planet is accreted onto 
the star (astrocentric distance less than 0.005AU), or 3) a planet is ejected from the system 
(astrocentric distance exceeds lOOAU). We manually verified that for the vast majority of 
systems not identified as unstable, the final orbits are qualitatively similar to the initial 
conditions. 



2.1.2. HD 187123 



Butler et al.l (119981 ) announced an Msinz = O.SMj^p planet in a 3-dav orbit around 
HP 18 7123, a close solar analog (M* = I.IM©, [Fe/H] = +0.1, %s = 5810). Iwright et al" 



(120071 ) announced the existence of an outer companion with orbital period > 10 yr. At that 
time, the radial velocity history was too incomplete to fully determine the orbit, which had 
not yet closed. I Wright et al.l (j2007f ) were nonetheless able to constrain the minimum mass 
of this 'c' component to be planetary (l.SMj^p < Msinz < lOMj^p). Since then we have 
obtained additional observations of HD 187123, and have found that the orbit has closed 
just recently. 

Assuming that there is no linear trend or detectable third planet in the system, the data 
constrain HD 187123c to have P = 10.4±1.2yr, e = 0.25±0.03, and M sin i = 2.0±0.3Mj^p. 
We present the newly determined orbital parameters in Table [Hand the latest radial velocities 
in Table [2] and Figure |H Figure [5] shows in P — M sini space, demonstrating that, despite 
our having observed somewhat less than one full orbit, the minimum mass and orbital period 
of the planet are constrained to ~ 20%. 

We have checked this orbit for long-term stability in the same manner as for HD 183263 
and find that the range of solutions shown in Figure [5] are well within the stable regime. 



2.1.3. HD 11964 



HD 11964 (= GJ 81. lA) is a metal-rich ([Fe/H] = +0.12) slightly evolved ( 2 mag, above 
the na ain sequence) G star with a nearby (sep ~ 30"), K dwarf companion. iButler et al. 



( 120061 ) announced the planet HD 119646, a Jovian (Msini = 0.6Mj^p) plan et in a 5.5 yr. 



circula r orbit, and noted a weak, residual trend in the velocities. An analysis by lWright et al. 
( 120071 ) showed that a trend was probably not the proper interpretation of the residuals, and 
that they we r e con sistent with a low amplitude {K =5.6 ms~^), 38 d signal (FAP < 2%). 
Wright et al.l (120071 ) cautioned that the low amplitude of this prospective 38 d planet meant 
that it would require more observations for confirmation, especially given the higher levels 
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Fig. 4. — RV curves for HD 187123. The data are from Keck Observatory and show the 
inner planet (top) with P = 3.1 d and Msini = O.SMj^p and the outer planet (bottom) 
with P = 10.7 yr and Msini = 2 Mj^p. 
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Fig. 5. — Contours of in Pc — rric sinic spac e of best-fit orbits to the RV data of HD 187123 
(Fig.H]), with in grayscale (cf. Figure 15 of IWright et al.l (120071 )). The solid contours mark 
the levels where increases by 1, 4 and 9 from the minimum. Although the orbit is still 
imprecise, the orbital period and minimum mass are now constrained to ~ 20%. 
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Fig. 6. — RV curves for HD 11964. The data are from Keck Observatory and show the inner 
planet (bottom) with P — 38.9 d and Msini = SSMq, and the outer planet (top) with 
P = 5.5 yr and Msini = 0.6 Mj^p. 



of jitter seen in subgiant stars (IWrightl l2005l : I Johnson et al.l l2007bl ). iGregoryl (120071 ) also 
noted the 38 d period in a Bayesian periodgram of the pubhshed velocities (as well as a 360 
d signal which is not apparent in the rereduced data presented here). 

We have now obtained 24 additional observations of the star, and the 38 d signal has 
strengthened (FAP < 1%), allowing us to confirm a low-mass (Msini = 23M0) planet, HD 
11964c, in a 37.9 d orbit. The low amplitude of the signal makes estimation of e and uj 
difficult, but data are inconsistent with a circular orbit, and favor e = 0.3. 

Private communications from our group regarding the 38 d signal, which has been 
apparent, but not convincing, since 20 05, has led to sorae cori fusion in the literature regarding 
the nomenclature of these planets (e.g. iRaghavan et al.ll2006l ) . Here, we follow the convention 
that planet components are ordered by the date of a formal or a public announcement of 
their existence. We report the orbital parameters of the two-planet fit in Tabled], the radial 
velocities in Table [21 and we show the radial velocity curves in Fig. [61 

We have checked this orbit for long-term stability with the same n-body code as for HD 
11964, and find that it is well within the stable regime. 



2.14. HD 217107 



Fischer et al. 



and Fischer et al. 



velocities. IVogt et al. 



1999|) reported a 7.1 d, Msini = 1.4Mj^p planet orbiting HD 217107, 



(l200ll) desc ribed a linear trend superimposed on the Lick and Keck radial 



(I2OO5I ) updated the orbital fit, finding significant curvature in the 



residuals to the inner planet fit, and estimated its orbital period, though poorly constrained, 
to be 8.5 yr, with Msini ~ 2Mj^p. 

Since then, data collected at Lick and Keck continue to map out the orbit of the outer 
planet. Today, we can constrain the minimum mass and period to within ~ 10% under the 
assumption that there are no additional planets in the system. 

We present the newly determined orbital parameters in Table [H and the latest radial 
velocities in Table [2l and Figure [3 Figure [H shows in P — Msini space, demonstrating 
that, despite our having observed somewhat less than one full orbit, the minimum mass and 
orbital period of the planet are constrained to ~ 10%, at P ~ 11.7 yr and M sini ~ 2.6Mj^p. 

We have checked this orbit for long-term stability in the same manner as for HD 183263, 
and find that the range of solutions shown in Figure [8] are well within the stable regime. 
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Fig. 7. — RV curves for HD 217107. The data are from Keck Observatory and Lick Observa- 
tory, and show the inner planet (top) with P = 7.1 d and Msini = 1.4Mj^p and the outer 
planet (bottom) with P — 11.6 yr and Msini = 2.6 Mj^p. 




Fig. 8. — Contours of in Pc — rric sinic spac e of best-fit orbits to the RV data of HD 217107 
(Fig. [7]), with in grayscale (cf. Figure 15 of IWright et al.l (120071 )). The sohd contours mark 
the levels where increases by 1, 4, and 9 from the minimum. Although the orbit is still 
imprecise, the orbital period and minimum mass are now constrained to ~ 10%. 



2.1.5. J^l UMa 



Butler &: Marcyl (119961 ) announced the existence of a 1090-day planet orbiting 47 UMa 



from data co l lected at Lick Observatory. After collecting an additional 6 years of data, 



Fischer et al.l (j2002al ) a nnounced th e exist ence of a second, 0.46 Mj^p long peri o d corn panion 



m a ~ 2600-day orbit. iNaef et al.l (120041 ) and IWittenmver. Endl. fc CochranI (120071 ) . using 
data from ELODIE and McDonald Observatory, respectively, have questioned the existence 
of 47 UMa c. Neither of the latter data sets, however, have both the precision and the 
duration to rule out the outer planet. The parameters quoted here are the literature values. 



2.2. The Current Sample 



We consider here the 28 known multiple-planet systems among the 205 known, normal 
exoplan et host stars withi n 200 pc. This is the sample of the Catalog of Nearby Exoplanets 
(ONE, butler et al.l [2006llN. The CNE employs a liberal upper mass limit in its definition 
of an exoplanet (any companion with Msinz < 24Mj^p) but restricts itself to systems with 
high-quality radial velocity detections around the bright stars most amenable to confirmation 
and follow "-up. This distance cutoff excludes the multiple planets orbiting the pulsar PSR 
1257 +12 (|Konacki et all WM) and the Jupiter-Saturn analogs orbiting OGLE-2006-BLG- 



109 (iGaudi et al 



20081 ) detected by microlensing. We also exclude several speculative claims 



of additional planets around known exoplanet host stars, systems where second planets have 
very poorly constrained orbits, and two announced multiplanet systems (HD 47186 and HD 
181433 by Bouchy et al. 2009) for which orbital parameters were not available at submission 
time. 

We provide up-to-date fits with radial velocities as recent as 2008 June for those mul- 
tiplanet systems with no significant planet-planet inte ractions (see [12.41) and for which we 
have Lick and Keck data from our planet search (see iButler et al.l l2006l . for details). We 
have employed a fitting algorithm which exploits linear parameters in the Kepler problem 
(Wright & Howard, 2008, ApJ submitted) to efficiently search the high-dimensional space 
associated with multiple-planet systems. We have also updated some of our radial velocity 
data reduction procedures at Lick and Keck Observatories, including a small correction to 
our calculation of telescopic barycentric motion. The radial velocities presented here are 
thus more accurate and precise than our previously published velocities for these systems. 
The resulting best-fit orbital parameters and uncertainties supersede previously published 



^Available on the World Wide Web at http: //exoplanets. org 



parameters. 

Table [H contains measured propert ies and derived quan tities for these 28 systems. Prop- 
erties of the host stars can be found in iButler et al.l (120061 ). 



2.3. The Multiplicity Rate 



In additio n to these 28 system s, 36 single-planet systems are best fit with the addition of 
a linear trend (IButler et al.ll2006h . If we conservatively exclude the eight such cases in which 
the host star has a known stellar companion (determined from a survey of the literature 
including lEggenberger et al.l 120071 ) . then we are left with 28 apparently single-star systems 
with an outer, potentially planetary companionjfl Of the 205 known nearby planetary systems 
then, 14% have multiple confirmed planets and another 14% show significant evidence of 
being multiple, meaning the tr ue planet mu l tiplici ty rate may be 28% or higher. This is 



Wrig 



it et 



consistent with the es timate of 
the value of ~ 50% in IPischer et al.l (l200lf ). 



al.l (I2OO7I ) and somewhat more conservative than 



2.4. Kinematic vs. Dynamical Fits 

The radial velocity signature of multiple-planets is significantly more complex than that 
of a single-planet system, and fitting such curves to observed radial velocity data requires 
care. Each planet has five spectroscopic orbital parameters, most of which are neither 
orthogonal nor linear, so finding a global minimum in space becomes significantly more 
difficult as the number of planets grows. Efficient algorithms are necessary to conduct a 
thorough search (Wright & Howard 2008, ApJ submitted). 

Planets in mean-motion resonances can be particularly difficult to identify from radial 
velocity curves because of degeneracies among the best-fit orbital parameters. Care must 
also be taken not to confuse weak signals resulting from al iasing of orbital pe r iods with the 



obser ving window function with genuine planet detections (jPischer et al.ll2008l ; iTinney et al. 



20061). 



Most importantly, interactions between planets may require consideration. In many 
cases these planet-planet interactions are sufficiently weak that they can be ignored, and the 



•^Our radial velocity analysis is sensitive to the presence of a second set of spectral lines, and we esti- 
mate that we can rule such a binary companion in a close orbit down to ~ 0.1 Mq in most cases. Some 
contamination from as- yet undetected binary companions may still remain in the sample, however. 



resulting radial velocity signal is simply the linear superposition of multiple Keplerian radial 
velocity curves (a "Keplerian" or "kinematic" fit). 

When these interactions are important, however, both short term and long term numer- 
ical n-body integrations of the physical system must be performed. In th e short term, these 



interactions can ca use detectable variations in the orbits of the planets (IRivera et al.l 12005 



Fischer et al.ll2008l ). In these cases a set of constant Keplerian orbital elements is insufficient 
to model the observed radial velocities, and a proper fit must be driven by an ra-body code 
(a "Newtonian", or "dynamical", fit). 

In some cases even good Newtonian fits to the data may yield orbital parameters for 
planets which, while stable for the duration of the observations, are not stable on timescales 
comparable to the age of the planetary system. Thus, long-term stability is an additional 
constraint that multiplanet fits must satisfy. 

The new fits listed in Table [1] are all Keplerian fits which have been confirmed stable 
with the n-body code described in §2.1.2[ The compiled literature fits are a mixture of 
Keplerian and Newtonian fits, and in some cases a more sophisticated Newtonian fit may be 
superior to t he published one. For instance, the orbits of the planets in the HD 74156 system 
reported by iBean et al.l (120081 ) a re apparently unstable on 10^ yr timesc ales, but an orbit 
with similar parameters is stable (IBarnes. Gozdziewski. fc Raymond! l2008l ). Our conclusions 
on the statistical properties of planets are insensitive to these details, but detailed work 
on planet-planet interactions and resonant dynamics should employ the published RV data 
directly, rather than the orbital elements presented here. 



3. Statistics of Multiple-Planet Systems 

3.1. Semimajor Axis Distributions 

Figure [9] shows the distribution of semimajor axe^for multiplanet and apparently single- 
planet systems. The semimajor axes of multiplanet systems appear to show some significant 
departures from the single-planet systems. Most strikingly, the pile up of Hot- Jupiter planets 
between 0.03 AU and 0.07 AU and the jump at lAU observed in the distribution for single- 
planet systems are both absent from the corresponding distribution for multiplanet systems, 
which appears rather uniform. 



"^Following iButler et al.l (|2006f ). we calculate semimajor axes from the measured orbital periods using 
Newton's version of Kepler's Third Law; the major source of uncertainty in a is usually the host star's mass. 



It thus appears that planets in multiple systems are not consistent with having been 
drawn randomly from the population of apparently single-planet systems. To provide a 
numerical measure of the difference between the two distributions, we created the cumu- 
lative distribution function (CDF) of the semimajor axis distribution (with an upper limit 
of ISMj^p) and applied a Monte Carlo Kolmogorov-Smirnov (K-S) test for a difference in 
distribution. The K-S test uses the two sided D-statistic, defined as the maximum difference 
between the two CDF curves. The Appendix contains a complete description of our Monte 
Carlo K-S procedure. 

Figure [TO] shows the CDF for these two distributions. A K-S test of the two distributions 
shows that we can reject the null hypothesis that the two samples were drawn from the same 
distribution with 90% confidence. 

This is somewhat surprising, since presumably many apparently single-planet systems 
have as- yet undetected low- mass planets and are, in fact, multiple systems. Presumably, 
then the difference is due to the presence of multiple giant planets within 5 AU. 

In particular, close-in planets (a < 0.07 AU) are not as common in multiplanet systems 
as in apparently single-planet systems. Among the 88 single planets having a < 1 AU, 38 
of them (43%) have a < 0.07 AU. A similar rate for the 39 planets in multiple systems with 
a < 1 AU would yield 15 planets within 0.07 AU, but only seven are known. 

The 1 AU jump in the single-planet sample may relate to the location of the "ice line" 
beyond which ices can participate in planet formation in the protoplanetary nebula. If this 
peak therefore represents planets which have not experienced significant migration, then 
the lack of such a peak in the multiple-planet distribution suggests that planets in multiple 
planet systems generally undergo more significant migration. 

Alternatively, the lAU jump may be indicative of a pile- up of planets at the typical 
orbital distance within which migration beco mes inefficient du e to either the presence of 



the ice line or a "dead zone" near lAU (e.g. Ilda fc Lid |2008| ). For instance, in a model 
where planet migration is driven by planetesimal scattering, the disk surface density of 
planetesimals is large enough to drive migration beyond the ice line, but not inside it. 



3.2. Msinz 

The distribution of minimum masses (Msini) of planets in multiple planet systems 
are shown in Figure [11], along with that for single planets, multiplanet systems exhibit an 
apparent overabundance of planets with Msini between 0.01 to 0.2 Mt„„, but this may 
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Fig. 9. — Distribution of semimajor axes of exoplanets for multiple-planet systems (solid) 
and apparently single systems (dashed). Note the enhanced frequency of hot jupiters and 
the jump in abundance beyond 1 AU in the single-planet systems. 



be amplified by a selection effect. When we find a planet around a star we tend to observe 
that star more frequently — making it more likely that we will find another planet that was 
not detectable beforehand. This appears to explain the detection of very low mass planets 
around 55 Cnc, GJ 876, and /i Ara. In these systems, more massive planets were known in 
advance, and an especially large number of observations were made to refine their orbits. 
The lowest-mass planets were found in the course of these detailed observations. 

The M sin i CDF of single-planet systems (Figure [T2l) is relatively featureless with the 
logarithmic Msini axis, but the CDF of multiplanet systems deviates markedly from its 
single-planet counterpart. We calculate = 1.36 and D~ = —0.089, so we find p{D~^) = 
1.8%, so we can reject the null hypothesis that our samples are from a common distribution 
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Fig. 10. — CDF of semimajor axis for known multiplanet systems (solid) and apparently 
single systems (dashed). 

with > 95% confidence. 

3.3. Eccentricity 

The distribution of eccentricities for single and multiplanet systems are shown in Fig- 
ure [131 Note that we have excluded planets with a < 0.1 AU from consideration here to 
remove the effects of tidal circularization on the analysis. Of the planets selected, both single 
and multiplanet systems exhibit a wide range of eccentricities from 0.0 to 0.8. The mean 
and standard deviation for single planets are 0.30 and 0.24, respectively, and 0.22 and 0.17 
for multiplanet systems. Also, 11 single planets have eccentricities above 0.7 (7% of this 
sample), but none of the multiplanet systems have an eccentricity above 0.7. 
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Fig. 11. — Distribution of Msinz of exoplanets, with an upper limit of 13Mj^p for known 
multiple-planet systems (solid) and apparently single systems (dashed) 



To provide a numerical measure of the significance of the difference between the two 
distributions, we computed the cumulative distribution function of the eccentricities and 
applied a Monte Carlo K-S test for difference of distribution (Figure [T^ . The eccentricity 
CDF for multiplanet systems is greater than the eccentricity CDF for single-planet systems 
at any given eccentricity, suggesting that multiplanets systems have systematically lower 
eccentricity. The K-S statistics are = 1.34 and D~ = 0.265, and p{D^) =< 1%, so we 
reject the null hypothesis that our samples are from a common distribution with over 99% 
confidence. Thus it appears that the known multiplanet systems have systematically lower 
eccentricities. Sele ction effects in e c centri city do not significantly affect the delectability of 
planets for e < 0.7 lCumming et al.l (120081 ). 



It is surprising that multiplanet systems have lower orbital eccentricities, as mutual 




Fig. 12. — CDF of Msinzfor known multiplanet systems (solid) and apparently single sys- 
tems (dashed) Note the enhanced frequency of known low mass planets in multiplanet sys- 
tems. 



interactions between giant planets might be expected to excite eccentricities. The lack of 
very high eccentricities in multiplanet systems may be partially explained by the additional 
constraint in multiplanet systems of orbital stability, which favors low-eccentricity orbits. 
Conversely, some single-planet systems may exhibit high eccentricities as a result of a series 
of ejections of former rnember s from the system. Both factors can be at play simultaneously: 



Ford. Lystad. fc Rasid (120051 ) explain the observed eccentricities of the planets in the v 



Andromedae system as the end result of the ejection of a hypothesized fourth planet from 
the system. 
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Fig. 13. — Distribution of eccentricities of exoplanets for known multiple-planet systems 
(solid) and apparently single-planet systems (dashed). Note the high eccentricity orbits, 
e > 0.6 occur predominantly in single planets. 



3.4. Metallicity 



Fischer fc Valentil (120051 ). using uniformly calculated metallicities from lValenti fc Fischer 



(120051 ). found that the 14 multiple-planet systems then known had a somewhat higher aver- 
age metallicity than single-planet systems (+0.18 versus +0.14), suggesting that metallicity 
trac es multiplicity in planets even mo re strongly than it traces single-planet occurrence (see 
also ISantos. Israelian. fc Mayorl l200ll ) . They suggested that further discoveries of multi- 
planet systems could confirm this trend. Using the metallicities compileclf] in the CNE, we 
find that single and multiplanet systems have mean [Fe/H] values of +0.10 and +0.10, re- 



■^In most cases, the ultimate origin of these metahicites is 
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Fig. 14. — CDF of eccentricity for multiplanet systems (solid) and apparently single-planet 
systems (dashed). The tidally circularized hot jupiters have been removed. Note that the 
highest eccentricities, e > 0.6 occur predominantly in the single-planet systems. 



spectively, although the median values are +0.15 and +0.18, still showing some evidence of 
the disparity. 

If we include systems showing long-term RV trends (not including known binaries) 
among the multiple-planet systems, however, the difference becomes slightly stronger. Such 
systems have a mean [Fe/H] value of +0.20, bringing the average for apparent multiple 
systems overall up to [Fe/H] = +0.15. 

We plot the two distributions in Fig. [151 To test whether these two distributions differ 
significantly, we have performed a K-S test, as shown in Figure dni A K-S test rejects, 
with 97% confidence, the null hypothesis that the apparently multiplanet systems (including 
systems with trends) have metallicities drawn from the same distribution as the single-planet 



systems. It appears that metallicity traces not only planet occurrence rate, but multiplicity 
among planet-bearing stars, as well. 
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Fig. 15. — Distribution of [Fe/H] for exoplanet bearing stars harboring likely multiple-planet 
systems (including single-planet systems with long-term RV trends; solid line) and apparently 
single systems (dashed). The median [Fe/H] for known multiplanet systems, +0.18, is higher 
than that for the single-planet systems, +0.14. The multiplanet system HD 155358 has the 
lowest [Fe/H] of any system, with [Fe/H] = —0.68. 



3.5. Stellar Mass 

The distributions of stellar mass for stars hosting single planets and stars hosting multi- 
planet systems are shown in Figure [T71 multiplanet systems have a mean stellar mass of 1 . 1 
Mq and single- planet systems have a mean mass of 1.13 Mq. K-S tests of the histogram in 
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Fig. 16. — CDF of [Fe/H] for known multiple-planet systems (solid) and apparently sin- 
gle systems (dashed). The CDF shows that the metallicity of the multiplanet systems is 
consistently higher than that of the single-planet systems. 

Figure [T7| and of the corresponding CDF show no significant difference between the stellar 
masses of single and multiplanet systems. Unlike metallicity, stellar mass does not seem to 
strongly trace multiplicity. 

3.6. Multiplicity vs. Stellar Mass 

Although there is no strong stellar-mass-multiplicity relation, there may be an emerging 
trend regarding the M dwarfs. Among the eight known M-dwarf exoplanet hosts, two are 
well characterized multiple-planet systems (GJ 876 and GJ 581) and two others (GJ 317, and 
GJ 849) show a trend in the single-planet velocity residuals. If we attribute these trends to 
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Fig. 17. — Mass distribution of exoplanet-bearing stars for known multiple planet systems 
(solid) and apparently single systems (dashed). There is no significant difference between 
the stellar masses of the single and multiplanet systems. 



planets and not undetected stellar or brown dwarf companions, then the total, true multiple- 
planet rate is at least 50%, which is higher than the 28% similarly calculated for the entire 
CNE sample in §2.31 It appears that planets around M dwarfs may be found preferentially 
in multiplanet systems rather than singly. This trend, if it is not simply Poisson noise and 
should hold up as more M-dwarf systems are discovered, is especially surprising since M 



dwarfs have a lower-than -average planet occurrence rate (lEndl et al.l l2006l : iJohnson et al. 



2007al : iMavor et al.l 120081 ). 



There may be subtle observational selection effects at work here, however. If M dwarfs 
in general have lower mass planets than F-K stars (and thus require more observations 
before publication), then we may simply be seeing the already-documented increase in planet 



occurrence rate amongst low-mass planets. The ongoing RV M dwarf surveys will improve 
the statistics of these systems, which should help illuminate if the effect is due to selection 
effects, small numbers, or astrophysics. 



3.7. Eccentricity vs. Msini 

Consider the plot of eccentricity vs. Msini of planets with a > 0.1 AU, including both 
multi and single-planet systems (Figure [T8|) . We have excluded planets within 0.1 AU to 
remove the effects of tidal circularization on our analysis. Planets with minimum mass below 
1.0 Mj^p have a mean eccentricity of 0.19, while planets above this threshold have a mean 
eccentricity of 0.34. 

Figure [19] shows the strong dichotomy between the eccentricity distributions of super- 
and sub-Jupiters: the eccentricity of sub-Jupiters peaks at e < 0.1, while the eccentricity 
of super- Jupiters is distributed broadly from 0.0 < e < 0.6. Figure [20] shows that a K-S 
statistic bears out this difference: for these two populations = 2.06 and D~ = —.034, 
yielding p{D~^) < 0.1%, so we reject the null hypothesis that the samples are drawn from 
the same distribution, and conclude that the apparent difference is not a chance result of 
small-number statistics. 

This correlation between eccentricity and minimum planet mass provides a valuable clue 
about the origins of eccentricities. One possibility is that conditions that encourage the for- 
mation of planets with Msini > l.OMj^p may contribute to greater eccentricity pumping as 
well. If so, then these would need to be strong effects as more massive planets are rela tively 



more difficult to perturb. Also possibly relevant is the work of iGoldreich &: Saril (l2003l ). who 
describe a mechanism for modest eccentricity pumping of a planet through interactions with 
a protoplanetary disk that should be more efficient for more massive planets. This mecha- 
nism could create a "seed" eccentricity preferentially in massive planets, which would grow 
through planet-planet interactions. Alternatively, it is possible that significant eccentricities 
are the norm for all systems shortly a fter the diss ipation of the protoplanetary disk triggers 



strong planet-planet interactions (e.g. lFord]l2006l ). and that circularization through dynam 



ical friction with planetesimals is more efficient for planets with M < l.OMj^p. Another 
possibility, that this is the signature of eccentricity pumping by planet-planet scattering, is 
discussed in §3.10[ 
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Fig. 18. — Plot of eccentricity vs. M sin i for planets with a > 0.1 AU, to avoid contamination 
from tidal circularization. Filled circles are single-planet systems, and open circles represent 
multiplanet systems. There is a slight sense of increase in the upper envelope in the range 
0.05 - 1.0 Mj^p 



3.8. Eccentricity vs. Metallicity 



We plot eccentricity vs. metallicity ([Fe/H]) for planets with a > 0.1 AU, including both 
multi and single planet systems (Figure [2T1) . There appears to be no significant correlation 
between [Fe/H] and eccentricities. While the systems with [Fe/H] < —0.2 appear to have 
lower eccentricities than those with [Fe/H] > —0.2, this may be due to the small number of 
such systems. 
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Fig. 19. — Distribution of eccentricities of exoplanets with Msinz < l.OMj^p (solid) 
Msinz > l.OMj^p (dashed). The tidally circularized hot jupiters have been removed. Note 
that the eccentricity of planets of minimum mass < l.OMj^p peaks at eccentricity < 0.2, 
while the eccentricities e of planets of minimum mass > l.OMj^p are distributed broadly 
from 0.0 < e < 0.6. 



3.9. Eccentricity vs. Semimajor Axis 



Figure [22] shows the eccentricity vs. semimajor axis for all planetary systems. There 
is a clear pauci ty of planets with small a and large e, as e xpected from the effects of tidal 
circularization (IRasio et al.lll996l : iFord. Rasio. fc Sillslll999l ). 



Planets from 0.5 AU to about 3 AU have the widest range of eccentricities, and beyond 3 
AU there may be a paucity of planets with e >0.6. This could easily be due to observational 
biases: planets beyond 3 AU have such long orbital periods that only a small number of 
orbits have necessarily been observed. In cases with fewer than two complete orbits, the 
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Fig. 20. — CDF of eccentricity for planets of minimum mass < l.OMj^p (solid line) and 
minimum mass > l.OMj^p (dashed line) The tidally circularized hot jupiters have been 
removed. Note that the CDF for the hght planets rises steeply from eccentricity 0.0 to 0.2, 
while the CDF for the heavy planets rises with a nearly uniform slope from eccentricity 0.0 
to 0.6. 



radial velocity signature of any additional exterior planets with orbital periods longer than 
the span of the observations can sometimes be absorbed into the eccentricity term of the 
orbital solution of the inner planet. Further, since an e > 0.6 planet spends a small fraction 
of its orbit near periastron, where its velocity signal is largest, such a planet may not reveal 
itself until it has completed nearly an entire orbit. The discovery of more long-period planets 
and the observation of more complete orbits will help reveal the true eccentricity distribution 
of planets with a > 3 AU. 
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Fig. 21. — Plot of eccentricity vs. [Fe/H]for planets with a > 0.1 AU, to remove bias caused 
by tidal circularization. Filled circles are single-planet systems, and open circles represent 
multiplanet systems. No correlation is apparent. 



3.10. Ratio of Escape Speeds vs. Eccentricity 



Some theoretical models of the evolution of eccentricity through planet-planet scattering 
focus on the parameter 9, defined!^ as the ra tio of the escape speed from the planet to that 



of the planetary system (IFord fc Rasioll2008l ): 



GM 

Rr, 



^Some authors refer to this quantity as the Safronov number. 
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Fig. 22. — Plot of eccentricity vs. semimajor axis. Filled circles are single-planet systems, 
and open circles represent multiplanet systems. Note the increasing upper envelope from 
0.01 to 0.5 AU. We have drawn the curves representing periastron passage distances of 
a(l — e) = 5i?0 and a(l — e) = lOi?© with dashed lines for reference. 



where M is the mass of the planet, its radius, r its orbital distance, and is the mass 
of the host star. Because we do not have exact masses or radii for many of the exoplanets 
in our sample, and because the following analysis is rather insensitive to the exact values of 
those quantities, here we approximate 



^^sinz = 10 \\^\\^\[ (2) 




and crudely estimate exoplanetary radii from the assumption that exoplanets below the mass 
of Jupiter have similar mean densities: 



Rjnp) \ 1 Msini > Mjup 

When 6 > 1, a planet can efficiently eject bodies during close encounters, and when 6 < 1 
collisions are more frequent. Figure [231 shows the distribution of 6"^ si n i vs. eccentri c ity fo r 



single planets and those in multiple-planet systems. Consistent with iFord fc Rasid (l2008l ). 
Figure [251 shows that planets that scatter planets and planetesimals efficiently (i.e. those with 
6 » 1) have a wider range of eccentricities than inefficient scatterers. This is consistent 
with the hypothesis that planet-planet and planet-planetesimal scattering is a dominant 
mechanism for the excitation of exoplanetary eccentricites. 

Because 6"^ oc Msini, these results are consistent with §3.71 where we showed that the 
eccentricity distribution of planets with Msini < IMj^p is peaked at very small e: the 
eccentricity dichotomy between high- and low-mass planets may be related to the ability of 
high-mass planets to more efficiently scatter planets and planetesimals. 

If we divide the sample into those planets with 6 < 1 and those with 6 » 1, we can 
examine whether the eccentricity distributions of efficient scatterers in single and multiplanet 
systems differ from one another. Figures [241 fc [25] show that inefficient scaterers in single and 
multiplanet systems have very similar eccentricity distributions. Interestingly, Figures [26] 
& [27] show that efficient scatterers (those for which 6^ > 1) in multiplanet systems appear 
to have higher eccentricities than those in single-planet systems. A K-S test for these two 
populations gives = 1.36 and D~ = 0.11, yielding p{D~^) = 2%, showing that it is 
unlikely that these samples are drawn from the same distributions. This is despite the 
overall tendency for multiplanet systems to have lower eccentricities (see §3.31 ) 



3.11. Metallicity and Stellar Mass vs. Msini 

Since metallicity and stellar mass both correlate with the occurrence rate of planets, it 
is reasonable to check to see if either also correlates with minimum planet mass. In Figs. [28] 
& [29] we plot M sin i versus these two quantities. Neither figure shows a strong correlation. 
There is a dearth of low-mass planets found around high-mass stars, but this may simply be 
an observational artifact, as such planets would have a lower reflex amplitude and so be more 
difficult to detect. Likewise, the typical minimum mass of planets orbiting M-dwarf stars 
appears to be lower than that around solar mass stars, but this may simply be an artifact 
of the fact that these low-mass planets are more detectable around M dwarfs, and that M 
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Fig. 23. — Eccentricity vs. the ratio of tiie escape velocity from the planet to the escape 
velocity of the star [9'^). 
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Fig. 24. — Distribution of eccentricites for single- and multiple-planet systems with ^ < 1. 




Fig. 25. — CDF for eccentricity for 6 < 1 planets in single and multiplanet systems. The 
eccentricity distributions for these inefficient scatterers do not differ significantly. 
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Fig. 26. — Distribution of eccentricites for single and multiple-planet systems with 9 > 1. 




Fig. 27. — CDF for eccentricity for > 1 planets in single and multiplanet systems. The 
eccentricity distributions for these efficient scatterers differ significantly. It appears that 
efficient scatterers in multiplanet systems display larger eccentricities than those in single- 
planet systems. 



dwarfs appear to have a lower planet occurrence rate overall. 
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Fig. 28. — Plot of [Fe/H] vs. Msini. Filled circles are single-planet systems, and open circles 
represent multiplanet systems. 



3.12. Semimajor Axis vs. Stellar Mass 



We plot semimajor axis vs. stellar mass (Figure [30ll . Two features are readily apparent 
in this plot: a lack of close-in planets orbiting stars with M > 1.5Mq and a lack of long- 
period planets orbiting stars with M < 0.5Mq. Thus, the semimajor axes of giant planets 
correlate positi vely and sensitively wit h stellar mass. The first of these features has already 
been noted by iJohnson et al.l (l2007al ). who find that the effect is statistically significant. 
The lack of long-period planets around M dwarfs is puzzling, since significant numbers of 
M dwarfs have been a part of the major radial velocity planet searches since at least 1995, 
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Fig. 29. — Msini vs. stellar mass. Filled circles are single-planet systems, and open circles 
represent multiplanet systems. Apparently, minimum planet mass does not correlate well 
with stellar mass in the range of 0.3 - 1.9 Mq. 

sufficient to detect any massive long-period planets at orbital distances of a few AU. We are 

undertaking a more thorough study of the occurrence rate of long-period planets around all 
of our targets to confirm the reality of the apparent dearth of long-period planets around 
low-mass stars. Nonetheless, a correlation seems to be emerging, driven primarily by the 
lack of long-period planets orbiting M dwarfs. 

3.13. Mass Ratio 



For the multiplanet systems, we plotted the ratio of the minimum mass of the outermost 
planet to the minimum mass of the next outermost planet M sin i (outer) /M sin i (inner) as 
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Fig. 30. — Plot of a vs. stellar Mass. Filled circles are single-planet systems, and open circles 
are multiplanet systems. There is no strong correlation among planets orbiting stars with 
0.8 < Mq < 1.3, but planets around the lowest-mass stars have smaller semimajor axes. 



shown in Figure [SH The distribution peaks near 1 and appears somewhat skewed toward 
systems in which the outer planet is more massive. However, the dependence of orbital 
distance on reflex amplitude undoubtedly plays an important role, since low-mass planets are 
more easily detected closer to their parent star. While it is true that detected outer planets 
tend to be more massive, this may not reflect the actual distribution of planet masses in 
multiple-planet systems. 
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Fig. 31. — Distribution of Msini ratios computed as M sin i (outer) /M sin z (inner) for the 
outermost two planets in the 20 known multiplanet systems. 



4. Summary and Discussion 

A major outstanding question in planet formation remains: how common are planetary 
systems with architectures similar to the solar system? In the known planetary systems 
within 200 pc, the median Msini is 1.6 Mj^p, and the median orbital distance is 0.9 AU. 
Given the observational biases, these numbers are suggest that the solar system may yet 
prove to be typical. 

The primary way in which the solar system appears unusual is in the nearly circular 
orbits of its planets, since the median eccentricity of the detected exoplanets is 0.25 (excluding 
tidally circularized planets). Figure [181 however, shows that these high eccentricities are 
generally restricted to planets with Msini > 1Mt„„, and that nearly circular orbits are 



typical in planetary systems having no planets with M > IMj^^. 

We have also shown that there are four classes of planetary systems with distinct orbital 
distance distributions — 1) planetary systems around stars with M > 1.5Mq, 2) those with 
M < 0.5Mq, 3) other apparently sing le-ptaet .y.tems, and 4) multiple-planet =ystemfl 
This provides a challenge for planet formation and migration theories to reproduce these 
disparate distributions. 

The 1 AU jump in the orbital distance distribution (Fig. [9]) of planets orbiting single 
stars may correspond to the planet formation ice line, or at least constitute an observational 
constraint on its location. It also suggests the existence of a large, yet unobserved population 
of giant planets beyond 1 AU. The next decade of radial velocity planet detection will reveal 
the population of giant plants from 3 to 7 AU, and thus provide a census of giant planets 
which have not experienced significant inward migration. This provides theory with an 
opportunity not only to explain the above observations, but also to provide observational 
tests of those explanations by predicting the frequency of such in situ planets orbiting various 
classes of stars. 



A. Statistical Methods: K-S tests 

In this work, we used the Kolmogorov-Smirnov (K-S) statistical test, based on the two 
sided D-statistic, to measure the likelihood that any appa rent differenc es in distributions 



(for a given property) were in fact attributable to chance (IKnuthl 119971 ). Considering any 
independent variable, "x", we construct two CDFs of that variable, one for planets in single- 
planet systems, Fs{x), and another for planets in multiplanet systems, Fm{x). We define the 
D-statistic as follows: 

= max{Fm{x) — Fs{x)) 

= max{Fs{x) - F„(x)) (Al) 

This statistic measures the maximum difference between the two CDF graphs being com- 
pared, with measuring how high the multiplanet CDF rises above the single-planet CDF, 
and D~ measuring how high the single-planet CDF rises above the multiplanet CDF. A large 
implies that multiplanet systems systematically have lower values of x, and a large D~ 
implies that they have higher values. We quantify this using a Monte Carlo method to 
produce a K-S confidence value for this statistic. 



^There is, at present, only slight overlap between these classes: GJ 581 and GJ 876 are both multiple- 
planet systems and have M < 0.5Mq 



To establish confidence levels we count n planets in single-planet systems and m planets 
in multiplanet systems. The "D" statistic measures how frequently m planets drawn with 
replacement from the distribution of x for single-planet systems will have D~^imuiated ^^"^ 
^simulated large or larger than the D-statistics for the actual multiplanet systems. We 
simulated 100,000 multiplanet distributions using the distribution for single-planet systems, 
and the fraction of trials with D'^i^^iated gi'eater than the for multiplanet systems is the 
K-S test p-value (p) that the CDF for multiplanet systems hes above the CDF for single- 
planet systems - and similarly, the fraction of trials with Dgimuiated greater than the D~ for 
multiplanet systems is the K-S test p-value that the CDF for multiplanet systems lies below 
the CDF for single-planet systems. The K-S test does not assume any particular shape for 
the distributions, however, we have assumed that the counts follow Poisson statistics. 
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